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Abstract 



We study a two-phase modified Stefan problem modeling solid combustion and nonequilib- 
rium phase transition. The problem is known to exhibit a variety of non-trivial dynamical 
scenarios. We develop a priori estimates and establish well-posedness of the problem in 
weighted spaces of continuous functions. The estimates secure sufficient decay of solutions 
that allows for an analysis in Hilbert spaces. We demonstrate existence of compact attractors 
in the weighted spaces and prove that the attractor consists of sufficiently regular functions. 
This allows us to show that the Hausdorff dimension of the attractor is finite. 
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1 Introduction 



The subject of this paper is a study of dynamics of a nonequihbrium two-phase Stefan 
problem modehng condensed phase combustion and some phase transition processes. It 
was demonstrated numerically 0[ that the sharp interface model of the condensed phase 
combustion also known as Self-propagating High-temperature Synthesis (SHS) generates a 
remarkable variety of complex thermokinetic oscillations. In addition to its theoretical inter- 
est SHS finds technological applications as a method of synthesizing certain technologically 
advanced materials, see |2II, 120] and also |2ZI for a popular exposition. The process is 
characterized by highly exothermic reactions propagating through mixtures of fine elemental 
reactant powders, resulting in the synthesis of compounds. The dynamical scenarios exhib- 
ited by the model include a Hopf bifurcation, period doubling cascades leading to chaotic 
pulsations, a Shilnikov-Hopf bifurcation etc. These scenarios are well-known for the finite- 
dimensional dynamical systems and suggest a possibility that the essential dynamics of the 
free-interface problem may be finite-dimensional as well. Indeed, we have been able to prove 
[T^ IT^ that compactness and finite dimensionality of the attractor take place for a simpler 
one-phase problem. 

However, the methods of the papers dealing with the one-phase problem are not directly 
applicable to the sharp interface problem of the condensed phase combustion which is the 
subject of the present paper. The principal difficulty that arises here as compared to the 
one-phase problem is that the presence of the additional temperature field behind the 
propagating interface (in the product phase) creates an additional degree of freedom that 
is not easily controllable. This difficulty is overcome in the present paper; we show that 
Hausdorff dimension of the attractor is finite. The paper draws on the approach of our 
previous work J3] discussing compactness of the attractor; in addition to the dimension 
estimate, results presented in this paper clarify structure and regularity properties of the 
attractor. 

There is a substantial literature that treats analytical aspects of the initial-boundary 
value problem for different sharp-interface models with kinetics related to the problem (jl.H - 
II. 4|) below, see |22l El El ■ These works are concerned with basic issues of mostly local in 
time existence. We also note recent papers by Brauner et al. and Lorenzi, 01211^; which 
study weakly-nonlinear dynamical behavior of solutions of related problems. In particular 
they consider perturbations of traveling-wave initial data and investigate their instability and 
bifurcations. In contrast, the principal focus of the present paper is in strongly nonlinear 
asymptotic dynamics for a wide range of initial data and parametric regimes. 

The free-interface problem is formulated as follows: find s{t) and u{x,t) such that 



Ut = U: 



•XX 



7M, X 7^ s(t), t > 







where v{t) is the interface velocity, s{t) = v{T)dT is its position, u is the temperature. 
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and the derivatives and u~ are taken from right side and left side of the free interface 
respectively. The last term in the heat equation p.lj) is due to the heat losses into the 
medium surrounding the combustible or solidifying substance via Newton's cooling law with 
a non-dimensional coefficient 7 > 0. 

Dynamics of the physical system is determined by the feedback mechanism between 
the heat release due to the kinetics g(u\x=s{t)) and the heat dissipation by the medium. 
The second interface condition ()1.4p (the Stefan boundary condition) expresses the balance 
between the heat produced at the free boundary and its diffusion by the adjacent medium. 
As the problem describes, generally speaking, propagation of a phase transition front, the 
first interface condition is a manifestation of the nonequilibrium nature of the transition; 
its analog for the classical Stefan problem is just u\x=s(t) = 0. We should mention that in 
contrast with the nonequilibrium problem, the dynamics of the classical Stefan problem is 
relatively trivial. The surrounding matter is assumed to be at the temperature of the fresh 
combustible mixture at —00 (the original phase in the phase transition interpretation). By 
the same token the heat loss will reduce the temperature in the product phase to that of the 
medium. Thus the behavior of the solution at infinity should satisfy limx^±oou{x,t) = 0. 

In order to estimate the Hausdorff dimension of the attractor we need to develop some 
additional technical tools. We first develop a priori estimates and establish well-posedness of 
the problem in weighted spaces of continuous functions. These estimates, which constitute 
the analytical core of the paper, secure sufficient decay of solutions that allows us to carry 
out analysis in a Hilbert space. It should be noted that volume estimates, which are the 
basis for the Hausdorff dimension bound, require a Hilbert structure in the underlying space. 
Next we are able to extend our results on existence of compact attractors ^3] to the weighted 
spaces. After that we study the evolution on the attractor and prove that the semigroup 
on the attractor is onto and one-to-one: it yields that the attractor consists of sufficiently 
regular functions. As a consequence we are able to demonstrate that problem is well-posed 
and its attractor is precompact in a Hilbert space. 

This allows us to estimate the Hausdorff dimension of the attractor based on the tech- 
niques described for instance in [23] • We study evolution of the infinitesimal volume along 
the trajectories in the attractor and demonstrate that for sufficiently large m that is de- 
fined solely by the physical properties of the problem, the m-dimensional volume decays 
exponentially. After that we prove that the semigroup is uniformly differentiable which, 
combined with the estimate for the linearized evolution of the infinitesimal volume leads to 
the conclusion that the Hausdorff dimension of the attractor is finite. 

2 Properties of solutions: Previous results 

In this section we present some pertinent background information from i jj (certain statements 
are slightly modified and clarified). The following theorem summarizes existence results: 

Theorem 1 Suppose that the kinetic functions g satisfies the following assumptions: 
{Al) g{u) is a continuously differentiable, monotone decreasing, negative function on (0, 00) 
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with g{Qi) = —vq for some velocity —vq < 0; 
(y42) g{u) is suhlinear: \\mu->oogiu) /u = 0; 

and that the initial data uq{x) G C(— 00,00). Then there exists one and only one classical 
solution of the free interface problem M.l\) - \l.Ji\j . The solution is uniformly hounded for all 
t > 0. 

The proof is based on the reduction to an integral equation for the interface velocity 
v{t) = g f e^^* r G{s{t), t, e, O)uo{0d^ - [' G{s{t), t, s{t), r)e~"'^'~Mr)dr) , (2.5) 







which arises from the interface condition and the single layer potential representation for the 
solution operator: 



{Tuo){x, t) := u{x, t) = e"^* / G{x, t, ^, 0)^0(0^^^ - / Gi^, i, 4^), r)e~^(*"")f (r)rfr, 

J -00 Jo 

(2.6) 

where 

G(x,t,e,r) =exp{-|^}[47r(t-r)]-V2 ^2.7) 

is the heat kernel and s{t) = v{T)dT. 

In the sequel we replace the sublinearity condition {A2) by a stronger condition. We 
assume that g{u) is a monotonically decreasing differentiable function on [0, 00] with \g'\ < C 
and satisfying 

< g{u) < -Vq for some 1/°, Vq > 0. (2.8) 

These conditions are satisfied, for instance, for the standard Arrhenius kinetics where v = 
V^expi-A/{u-Uoo)). 

Under some additional conditions the following smoothness result holds [?]: 

Theorem 2 Let the initial data uq be twice differentiable in x < and x > with bounded 
derivatives and satisfy the matching condition: 

,(no(0)) = ^(0) - ^(0), (2.9) 
ox ox 

in addition, let the derivative of the kinetics function g' be Lipschitz continuous. Then the 
velocity v is differentiable. 

3 A priori estimates in weighted spaces 

For our purposes we need to establish certain a priori bounds on the solution in appropriate 
weighted spaces that are introduced next. Let tOa be the weight uja{x) = e"'^'; we define 

\f\a = sup{uJaix)\f{x)\), C„ = {/ G C(-oo, oo) : I / |„ < 00} 
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We will demonstrate that the global existence results can be can be extended to Ca (note 
the obvious imbedding, if /5 > a > then <Z Ca C C{— 00,00)). Similarly we define 
Hilbert versions of weighted spaces: 

11/ L = 11^"/ IIl2(-oo,oo) , Ha = {f ■■ 11/ IL < 

It is easy to see that if P < a then Ca C H[^, and for any f & Ca 

(poo \ 1/2 / poo ^2 \ 1/2 ^ 

Ij^lsH.) -[ljf.f^^) <-\n.^. (3.10) 

It is convenient to split the representation formula (j2.6p for the semigroup operator T 
into two operators: the contribution of the free boundary 



T^{t)u\x') = - e-'<^'-'^C{x\t,s{T),T)[v{T)]dT (3.11) 
Jo 

and that of the initial data 

/oo 
G{x\t,^,Q)u\i)di (3.12) 
-00 

In the sequel we will frequently encounter integrals of the error function type. To estimate 
them we employ the following simple result. 

Lemma 3 For a,b > 

A=exp(— 6a^), for a > l/y/b 



f ^ ^ 2^ , / ^exp -6a^ , /ora > 

Ja [ /or < a < 1/Vb 

Proof. If aVb > 1 then 

1*00 poo 1*00 

/ exp(—b'r]'^)dr] = —= / exp(— 77^)^77 < / r] exi)(—b'r]'^}dr] = — ^exp(— 6a^) 

Ja Vb JaVb Ja 2\/b 

On the other hand 



exp(— 6?7 )dr] < j exp(— 6// )dri = 

2v6 



3.1 Estimates for the solution: Contribution from initial data 

Proposition 4 For sufficiently small a (if a satisfies + aV^ — 7 < the contribution 
from the initial data in the Ca-norm decays exponentially in time: 

\u2i;t)\a < 2exp[(-7 + a2 + tt1/°)t] \Uo\a 
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Proof. For the contribution from the initial data, U2{.,t) = T2{t)uo, we have: 
\u2{.,t)\a = sup \uja{x - s{t)) \ / e'^^G{x,t, ^, 0)uo{^)d^\ 



•J —CO 

oo 



L J-oo ^aKK) 

<— =|«oUsupK(x-.(t))[/ e"«exp(-^— / e""« exp(-^— 

2\/tTT J-oo Jo 4t 

Each of the integrals should be estimated separately for x < and x > with the maximum 
of the estimates chosen for the estimate of the norm. At the same time it is easy to see that 
the two integrals can be transformed into each other through the change x —x, therefore 
it suffices to estimate only one of them and double the result. 
We proceed as follows 



^, e-"«exp(-^^-^)rfe = e"'*e-"M exp(-7]''}dr] 
For X > the above expression 

/oo 
exp{-r]'^}d7] < v^e"'*e-"^' 
{x — 2ta.) 



OO 



2v/t 



and therefore 



Uo\a sup uja{x — s{t))^/^Te°' *e = |mo|q supexp(— 7)f: + a^t — ax) exp(a[x — s{t)]) 

TT x>0 x>0 

2 „ , „taOn 



< |Mo|aexp(a — 'y + aV )t 

where we have used the bound —s(t) < VH. 
To estimate the integral 

POO 

e"'*e-"^ / exp(-7]Hdr] 

J (2ta-x) 

for X < 0, we apply Lemma El If the lower limit is larger than unity then 

.aH^-a. r ^2, / 1 „„„,„2. „^ (2^ - x)\ _ 1 x' 



e / exp(— r? }dri < - exp(a t — ax : ) = - exp(— — ) 

2y/t 

If the lower limit is no larger than unity then 

/ exp(— ?7^}(i?7 < exp(a^t — ax) 

J (2ta-x) 2 



2Vt 

a/vt {2ta — x)"^ x^ eJli x^, 

^exp(^^)exp(--) < ^exp(--. 
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Thus, for X < 0, we get 

— ^|Mo|aSupexp(Q;|a; - s{t)\ - — ) < 

2 2 

-e~'''*|Mo|a iiiax[ sup exp{— ax + as{t) — — }, sup exp{aa; — as{t) — — }] 

2 x<s{i) 4i s(t)<i;<0 4i 

= -e~'''*|ito|a niax[exp(a^i — ai>ot), exp(y°to)] 
In the above estimate we used the elementary inequahty: 

2 

-ax + as(t) - — = -(^= + \foLtf + a^t + asit) < 
a^i + as{t) < a^t — avot 

(note that s{t) < —Vot). 

Collecting the estimates for all the cases {x < 0, and x > 0) 

\u2{-^t)\a = l^iola^"^* max[- exp(a^t — afot), - exp{V^ta) ^ exp{a^ + aV^^)t] 
< ^ exp[(-7 + a^ + ay°)i] |wo|a < 2 exp[(-7 + a^ + ay°)i] IwoU 

Thus, for any a we have obtained an a priori estimate on the contribution from the initial 
data valid for all time. If a is sufficiently small, + aV'^ — 7 < 0, then the norm of the 
contribution is exponentially decaying. We also note that for a — > the estimate has a limit 
and takes the form 

g 

|M2(-,t)|o < -exp(-7t) |tio|o (3.13) 



3.2 Estimates for the solution: Contribution from the free inter- 
face 

Proposition 5 The Ca-norm of the contribution from the free interface is uniformly hounded 
for all time: 

|(Titxo)(.,t)|a<^7V7, 
provided a < agpace '■= iiiin(t'o/4, 7/(21^°)). 

Proof. To estimate the free-interface contribution to the solution T-i(t)uo behind the 
interface x > s{t) we split the interval of integration into two subsets: Xi = {r G [0,i] : 
s(r) < {s{t) + x)/2} and its complement X2 = {t e [0,t] : s(t) > {s{t) + x)/2}. 

mt)Uo\< J' G{x,t,s{T),T)e'^^'-^^ \v{T)\dT = j + j = /i + /2, 

XI X2 
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For the first integral we have 

exp[— (x — s(r))^ ../ J 



- 2^ J ^ ' ' ^ " 16(t-r)- 

XI 
XI 

^ (a;-s(t))/(2i,o) 



erf Jt^^ exp[-(x - < -^exp[-(x - s{t)f-^] 



2^ I V 2vo j ^ ^^'4^-2^ ^ '"A 

The following inequalities 



(t-r) - ' 2 ' (t-r) - ' 2 x- s{t) 

have been used to replace the exponent in the Gaussian kernel, which after that gave the 
exponential decay factor. We note that the estimate has a regular behavior at the limit 
7 — > giving the bound 

It is a manifestation of the fact that the heat loss is immaterial in the vicinity of the interface 
(cf. the next estimate which shows that the presence of heat losses is essential for decay at 
large distances from the interface). 

For the integral I2 we use Lemma El to obtain 

/• e-^(*--^ exp[- ^'^7'^^y ]|t;(r)Mr 

X2 

00 00 

J 2^717/ 0r J 

:exp(-7(x - s{t))/{2V^)), for 7(x - s{t))/{2V°) > 1 



< <; (3.14) 
-— , forO<7(x-s(t))/(2\/0)<l 
2^7 
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Thus for X > s{t) we obtain 



2^/7 



exp[-ix - sit)n] + 



\mt)uoix)\ < < 



_Klexp(-7(x-.(t))/(2n), for 7(x - s(t))/(2n > 1 
:exp[-{x-s{t)rf] + 



(3.15) 



2^7 



2^/7 

for < 7(x - s{t))/{2V^) < 1 



We note that upon multiphcation by the weight exp(a(x — s(t))) the right hand sides of the 
estimate ()3.15p are decaying exponentials, provided 



a < ctspace ■= min(t;o/4,7/(2\/°)), 
that attain their maximum j ^ at x — sit) = 0. Thus 

\Ti{t)uo{x)\ < F7v^ for X > s{t) 
Ahead of the interface x < s{t) we have: 

I / e-7(t-r)^^ v(T)dT\ 

. -,it-r) , [(x-s(t)) + (.(t)-s(r))]% Kr)Mr 



(3.16) 



< exp{- 



t;o|2; - s{t)\ \x - s(t)|% V 



2 yO rt 2 
}^ / exp{(-^-7)(t-r)} 
Jo 4 



At 



dr 



2v^ 



+ 47 



Thus, 



At 



Ti(t)Mo(a;)| < ^==exp{ — ^ ^ for x < s{t) 



(3.17) 



+ 47 ' 2 
Now it is easy to obtain the estimate for the norm: 



|(TiMo)(-,^)U < max( sup [ 



.Vo 



x</{t)V^^o +47 



exp(-(^-a)(.(t)-a:))],— )<— (3.1^ 
2 \/7 



By combining the estimates for the initial data contribution and that from the free 
interface we arrive at the following result: 
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Theorem 6 For < a < agpace 



//, in addition, 



I (Two) la < — + 2exp[(-7 + + aV'')t] \uo\^ 
V7 



a < aume := — {^/l + 4^i/iV^ - 1) 



(3.19) 



(here aume is the positive root of + aV'^ — ^ — 0) then the contribution from the initial 
data decays exponentially. 



For the future use we combine the bounds: 



(3.20) 



Remcirk 7 We note that for realistic problems j/V^ <^ 1 then aume ~ l/V^, thus both 
bounds are of the same order and in this case ccmin = Oigpace- 



3.3 Estimates for the derivative: Contribution from initial data 

We also need estimates for the spatial derivative of the solution. We start with the contri- 
bution from initial data. 

Proposition 8 For a < aume the Ca-norm of the derivative of the contribution from the 
initial data decays exponentially in time: 

\{T2U).{., t)\^ < \uo\a{^ + ?) exp[(a2 + aV"" - 7)^]. 



Proof. First we split the integral 



|(T2«),(.,t)|„ = sup|a;c,(a;-s(t)) | j e-^'G,{x,t,^,Q)uo{Od^\ 



sup < e ^*uja{x — s{t))\ 



G,{x,t,i,Q)uo{i)uJami\ 



-It 



2Vt7r 



/■oo 

c,s\xp{uja{x - s{t))[j e~"^ 
Jo 



d (x-0^ 

4^' 



+ 



J —00 



d_ 

dx 



exp(- 



{x - 0' 
it 



Each of the integrals should be estimated separately for a; < and a; > and the maximum 
of the estimates should be chosen for the estimate of the norm. At the same time it is easy 
to see that the two integrals can be transformed to each other through the change x — > —x, 
therefore it suffices to estimate only one of them and double the result. 



9 



We proceed as follows. For x > we integrate by parts 



-exp( 



a 



Jo 



-f 

J X 



-exp( 



d_ 



exp(- 



)di 



)di- / e-«exp(- 



X 



)de) + 2e— -exp(--) 



a 



V 



exp(— — )(i?7 — / e 



4t 



jT ^-aix+v) exp(-|-)d^^| + 2e-"^ - exp(- 



4t- 



< q; exp(Q;^f — ax 



) / exp(- 

J —X 



2+2 



It 

[f] + 2atf 



■)dr] + 2e' 



/oo 
exp(— ' ' — ^)c^?7 + Se" 

= Q;exp(Q;^t — ax)\/t7r + 2e" 

For X < we integrate by parts to obtain 



^0 



d^ 



x-if 



d^ 



e-"«^exp(-^^^)cie = -a 



Now we are ready to estimate the norm: 



-7t 



-"(^+«) exp(-|)de + exp(-^) < exp(-^) 



\{T2u)x{-it)\a < — ^|rto|a ms-xisup -j uja{x — s{t)){aVtn exp(a^t — ax) + 2e~"^^) > , 
sup lujaix - s{t)) exp(-^)) 1] 



For the term with x > in the above estimate, we have 



Q,sup < u)a{x — s{t))[aVt7r exp{a^t — ax) + 2e °^)] > 



< — 7=|mo|o exp(— Q;s(t))[2 + ay/tTT exp{a'^t)] 

I^Jt'K 

< \uo\a |-^exp[(Q;y° - 7)0] + | exp[(Q;' + aV"" - 7)i]| 
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For X < 0, we have 



e"^* , , r x^ 



\{T2u),:{.,t)\c, < —^\Uo\a sup <UJa{x - S{t)) exp{- — )) 



< — =|Mo|aiiiax[ sup exp{— ax + as{t) — — }, sup expjaa; — as{t) — — }] 
2v^7r x<s(t) 4i s(t)<a;<o 4i 

< — ^|uo|a max |exp[(Q;'^ — a^o — 7)^], exp[(Q;V^° — 7)^]} 
2v ^TT 

x^ 

In the above estimate we have replaced ax — — by its maximum aH. We collect the estimate 
for X > and x < to obtain 

IM.(.,t)|a < \uo\a |^exp[(aV° -7)t)] + (^ + ^)exp[(a2 + aV^o j 

2 a 

< |iio|a(^ + 77) exp[(a' + - 7)*]- (3-21) 

Thus, for any a we have obtained an a priori estimate on the contribution from the initial 
data valid for all time. If < a < ajirne then the norm of the contribution is exponentially 
decaying. ■ 

3.4 Estimates for the derivative: Contribution from the interface 

Proposition 9 The Ca-norm of the derivative of the contribution from the free interface is 
uniformly bounded for all time: 

\{T,uU.,t)U < M{vo,V'',a,^) 

provided that 

< a < min(— , ^ ^ - 



Proof. (Ahead of the interface). The estimate ahead of the front x < s{t) is treated 
as follows. We consider separately two cases: \s{t) — x| > 1 and \s{t) — x\ < 1. 
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For the case \s{t) — x\ > 1 



l-t t \ -{x-s(T)Y/i{t-T) 







^ I 1 {X - s{t)) ^_(^_s(r))y8(t-T) ^ ^-(x-s(T))V8(t-r) v{r)dT 

2{t-T){x-s{T)) ^ATi{t - r) 

^1 f 4/e ^^^p, [{x - sit)) + {s{t) ^ s{r))]\ ^ v{r)dr , 



< 



^s{t)-s{T) %{t-T) ^nit r) 

4yOg-^o|x-.(t)|/4 rt U - s(t))\ , 

- e0F^;o|s(i)-x|7o 8(i - r) ^ 

— 7= — i — / e ^' df] < -- 



e^/7^vo\s{t) — x\ Jq e vo\s{t) — x\ 

In the last estimate we used the following simple observations: ^e~^ < 1/e, for ^ = 
(x - s(t))^ 

-^^ — - — > 0, |s(r) - x\ > \s{t) - x\, \s(t) - x\ > \s{t) - s(t)\ > Vo\t - rl and sub- 

A{t - r) 

stitution 1] = \s(t) — x\(t — r)^^/^ to obtain the error function integral. 
For the less involved case \s{t) — x\ < 1 we proceed as follows 



2(i-r) ^4n{t - r) 
* - s{t)\ + \s{t) - s{t)\ e-(^-«W)V4(t-r) 



~ ' Jo 2(i - r) ^in{t - t) ^ ' ' 

A 7o 4 40F Jo y/(t - r) 



- 2 20F 7o ' - 2 ^ vo' 

Thus for x < s(t) we obtain 

|(riu).(x,t)|<<( ^r*^'^^^^^^^' (3.22) 

— (1 + — ) for < s{t) - x<l 
2 vq 
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Remark 10 The proof above shows that if \s{t) — s{t)\ > vo\t — t\ which holds if the basic 
assumption on the kinetics in i2.^) is satisfied, then the derivative ahead of the interface 
X < s(t) decays exponentially 

(J^-vo\x-s{t)\/4. 

\{Tiu)^ix,t)\ < 1^^^^^ _ ^1 ll^^llc[o,t] (3-23) 

The exponent — fo/4 can be improved to —vq/{2 + e) (at the price of increasing C ). 

Proof. (Behind the interface). Now we consider the domain behind the interface, 
x > s(t). We split the interval of integration into two subsets: Xi = {r G [0,t] : s(r) < 
+ x)/2} and its compliment X2 = {r E [0,t] : s(r) > (s(t) + x)/2}. 



muUx,t)\ < p-^^^^G{xXs{r),T^^^^^ I + 



'2, 



XI X2 



For the first integral we have 



1 f \X - s{t)\ {x - s{T)y -^(t-r)l I \\j I 

<Z1 / l^-^(^)l expi- ("-^(^^^^ e-^(^-W 
- 4v/?y (t_^)3/2e^Pi i6(t-r) '^^ 

XI 



{x-s{t))/{2vo) 



° ^(x-s(t))2«o 



where rj = {x — s{t))/y/t — r. 

To estimate the last integral we apply Lemma El to obtain 

I ^^exp(-(a; - s(t))i;o/8), for (a; - sit))2vo > 1 

[ 2^°, for < (x - s{t))2vo < 1 

For the estimate on the C^-norm, the function will be multiplied by the weight. We note 
that upon multiplication by exp(a(x — s{t))) the right hand sides of the top estimate is a 
decaying exponential, provided a < Vq/S, that attain its maximum 2V^/ a/tt at x — s{t) = 0, 
while the bottom term is bounded by 21^*^ exp(Q;/(fo/8)) < 2eV^. Therefore the contribution 
of Ji into the norm is bounded above, for example, by 6V^: Ii < 6V^. 
For the integral I2 

{x-s{T)f 

n.-.Mi-pi- 4(.-.) . 

J 2((-t) 2Vir(« - t) ' " 

X2 
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we use simple geometric considerations that show that in the domain X2 if ^(r) > x then 

\x-s{r)\ ^ \s{t)-s{r)\ ^ 



2{t-T) - 2{t-T) - 2 
while for s(t) < x 



\X — SiT] 



Therefore the estimate for the derivative in this case reduces to the estimate for the function 
itself fim|) and yields: 

r ^ exp(-7(a: - s(t))/(2V0)), for ^{x - s{t))/{2V') > 1 
2 I ^, for < 7(a: - s{t))/{2V^) < 1 

Similarly to the argument for Ji, one can see that if a < '~f/{2V^) then the contribution of 
I2 into the Cc-norm is bounded by y°e/(2^/7). 
Thus for X > s{t) we obtain 

sup \{T,u),{x,t) exp{a{x - < ^^6/(2^) + 6V' < V\2/^ + Q) (3.24) 

x>s{t) 

if a < mm{f, 

Similarly, if a < t>o/4 then by employing (|3.22p we see that 

sup \(Tiu)x:{x,t) exp{—a{x — s(t))\ < max[ ,— (IH )e°] 

x<s{t) ewo 2 vq 

Finally, for the norm we get 

muU;t)\a < V max[— , ^(1 + — )e", 2/^ + 6] := Mivo, «, 7) (3.25) 

eVo 2 Vq 

The estimate holds if 

< a < min(^, ^) := (3.26) 

(Recall that Tiu is the contribution from the free interface and therefore the absolute bound 
on its derivative is independent of the initial data). ■ 

Remark 11 The choice of the factor 1/2 for the split of the domain of integration into 
two parts above is rather arbitrary; by choosing a factor approaching 1 we can improve the 
exponent of decay in the subsequent result, the corresponding coefficient though will increase. 
Consequently, the value — in the definition of be improved to become 



4 + e 
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3.5 A priori estimates in and 

We collect the estimates ()3.2H) - ()3.25|) to obtain the following result: 
Theorem 12 

2 a 

|(Tmo).(., t)L < ^ + |wo|a(^ + -) exp[(aV^° \ o? - ^)t\ 

y/tn 2 

Remark 13 It is easy to demonstrate via integration by parts that if the initial data Uq G 
and satisfy the compatibility condition [{uo)x]x=o = fl'('Wo(0)) then the singularity in the 

above estimate will not take place. The estimate in this case reduces to the estimate for the 
solution through the derivative of the initial conditions. 

We note that a weaker result holds for Hilbert norms. 

Theorem 14 Let uq G Ca then the solution satisfies u{.,t) G C((0, oo), iJ^) where (3 < a 
and 

1 2 a 

II (T7io)(., t)||^,i < ^^{M + \uoU-= + -) exp[(ay° + - ^)t] 

\/(X — p y/tlT 2 

+ — + 2 \uo\a exp[(-7 + + aV^)t]} 

Proof. The proof is very simple since both Tuo{-,t) and (Tmo)x(-,^) G Hp by virtue of 
the imbedding estimate ()3.1()|1 . ■ 



4 Well-posedness 

In this section we prove that solutions of the free boundary problem ()l.ip - ()1.4|) depend 
continuously on initial data. This result is used in the sequel to demonstrate smoothness of 
the elements of the attractor. 

Theorem 15 In Ca, < a < agpace, the solutions of the problem depend on initial conditions 
uniformly continuously. More precisely, if {u{x,t), s(t)}, {u{x,t),s{t)}, < t < a are 
solutions with initial data vP,vP G Ca, where cr > depends only on the norm \vP\a and 
\vP\a, then for < t < a 

sup \V{t) -V{t)\ < clu"^ -u°\a, (4.27) 

0<t<(7 

\u{. - S{t),t) -U{. - S{t),t)\a < C\U° -U°\a (4.28) 

Remark 16 We state and prove continuous dependence on initial conditions only locally in 
time. The argument extending this result to any fixed time is based on the a priori estimates 
and follows closely the proof of global existence. 
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Remark 17 For simplicity of presentation we include the proof only for the uniform norm 
a = 0. The modifications for the case a > are rather routine but somewhat lengthy and 
follow along the similar lines. Everywhere in the proof below we use the notation \\.\\ = |.|o 

Proof. The proof consists of two parts. First we establish continuity of the interface 
velocity. We will establish first the estimate in ()4.27j) and then use it to derive ()4.28|1 . Let v 
and V be solutions of the integral equation in ()2.5|) with initial data and respectively. 



\v — v\ 



g\e-^' Gis{t),t, O)uo{0d^ - / t, s{t), T)e-'^'-^\{r)dr 

J-oo Jo 
poo pt 



+ L 



oo 
oo 







Gisit),t,^,0)uoiOd^- 

t, sir), T)e~^^''^^v{T)dr - I G{s{t),t, 5(r), r)e-^'^'-^^v{T)dr 

Jo 

To continue the estimate we employ a " coordinate descent" : 

/oo 
(G'(s(t),t,e,0)-G'(5(t),t,e,0))no(Orf^ 
-oo 

/oo 
G(5(t),t,e,0)(Mo(0-%(0)^^e 
-oo 

G{s{t),t,siT),T)e-^^'-^\v{T) -v{T))dT 

[G{s{t),t, s{t), r) - t, -sir), r)]e-^(*--)^(r)rfr 

:= Le-^^Di + Le"^*D2 + LD^ + LD^ 
To estimate the first summand 



+ L 



iG{sit),t, 0) - Gm,t, o))uo{m ■-- 



6Guod^ 



we note that by the mean value theorem, 



where 



dG s' — 

6G={s- ~s) — {s' - e, t, 0, 0) = {s- ~s)^G{s' - ^, t, 0, 0) 



s' = s'it,0, sit)<s'<:sit). 
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Thus 



|.' - i\G{s' - C, t, 0, 0) = (20F)-i |.' - e|r V2e-(^'-«^/^* 

= (2V^)-i(8i)V2l^_L|e-(^'-«'/«*2V2(2i)-V2e-(.'-0 

< 4t^/^ciG{s' - 2t, 0, 0) < Acit^/^G{s - ^, 2t, 0, 0) 
where ci = {2^/w)^^ m£Lx{xe^^^). Therefore 



1^1 



) J~oo 

1 nt POO 

<2c,t"^- \v-v\dT Gis-^,2t,0,0)\u\O\d^ 

I Jo J-oo 

sup \u°{^) \ \\v - v\\. 



f 

J — c 



Obviously, 



/oo 
G(5(t),t,e,0)|Mo(0-«o(OMe< ||wo 
■oo 



For the estimate of we replace the exponentials by 1 and integrate to obtain: 



Finally, for 



[\G{s{t),t, s{t),t) - G{§{t),t, ~s{T),T)]e-^^'-Mr)dT 
Jo 

< \\v\\ [ \AG\e-^^^-^^v{T)dTv{T) 
Jo 
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The estimations are quite elementary and are based on the mean value theorem. First we 
note that 

I AG| := |G(s(t), t, sir), r) - t, 5(r), r) | 

= |G(s(t) -s(r),t- r, 0,0) -G(5(t) -s(r),t- r, 0,0)1 

= \.,{t) - s{r) - (5(() - S(r))| \^(-s',t- r,0,0)| 
2{t - T) 

-.0.0)1 

<i||t.-{.|||,s'G(,s',(-T,0,0)| 

where t < r' < t and s' is between s(t) — s{t) and s{t) — s(r). Since 

|s'| < max{|5(t) - s(r)|, |s(t) - s(r)|} < Vo(i - r). 
and IGI < Co(t — r)"^/^ we get the estimate 

|AG| < CoVoWv - v\\{t - ry/^ 

Thus, for we get 

ID4I < C^Wv - v\\t^^'^ 
We collect the estimates for Di- to obtain 

11^^ - ^^11 < \\uo - Moll + lit; - v\\ + C4v - v\\t^/^ + C2t^'^\\v - v\\ 

From this inequality we see that for t < a, where a is small enough 

II ~ii ^ 11^0 -^^oll , ~ II 

- - l-C,tV2||nO||-C4t3/^-C,tV2 ^ 2||% - -oil 

We note that the value of a depends only on ||m°|| therefore this result can be interpreted as 
uniformly continuous dependence of v on the initial data on any ball \\u^\\ < R. 
Thus, the estimate ()4.27|1 has been demonstrated. Next we note that 

\u{x~ s{t),t) -u{x- s{t),t)\ (4.29) 
< \u{x — s(t),t) — u{x — s(t), t)\ + \u{x — s(t),t) — u{x — s(t),t)\ 

The second term in the inequality above can be estimated via the mean value theorem: 

\u{x — s{t),t) — u{x — s{t),t)\ = \ux{x — s',t)[s(t) — s(t)]\ 

< \ux{x — s' ,t)\t \\v — v\\ 
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where s' is the intermediate value between s{t) and s{t). 

The estimates for the first term in ()4.29|) are based on the maximum principle. After the 
free interfaces are determined, both u and u solve the heat equation off their respective 
boundaries. Their difference w = u — u satisfies the heat equation in each of the three 
domains 

X < s{t), s{t) < X < S(t), s{t) < X, T < (7 

here we assumed that s(r) < s(r) and that a is such that the inequality holds for all r < a. 
It is easy to estimate the boundary values of the difference 

\w{{s{t),t)\ = \u{s{t),t)-u{s{t),t)\ 

= \g~\vit))-g-\m-^As',tMt) - ~s{t)]\ 
< L\v{t) -v{t)\ + \u^{s',t)\t\\v -v\\ 

where again s' is the intermediate value between s{t) and s{t). Because of the a priori 
estimate on the derivative 

\Ux{s',t)\ <M + \U0\^ 

(see Theorem IT2|l we obtain 

\w{{s{t),t)\ < \\v-v\\ {L + Mt + cVi\\uo\\). 

A similar estimate holds for the other interface. Of course the initial data for w is equal to 
uq — uq. From the maximum principle for each of the three domains we obtain that 

\u{.,t) - u{.,t)\ < Ci \\v - v\\ + C2 \\uo - Moll 



5 Absorbing set and attractor 

In this section we use the estimates obtained above to establish existence of bounded absorb- 
ing sets and of the attractor which is compact in the weighted space of continuous functions. 
In order to establish compactness of the attractor we need to make use of the heat losses. It 
can be verified that most of the estimates and analytical properties of the solutions can be 
obtained without the heat losses. The presence of heat losses only improves the estimates. 
On the other hand the problem with the heat losses exhibits uniform exponential decay in 
time of the contribution of initial data which is utilized in the proof of compactness of the 
attractor. 

For our purposes we rephrase Theorem EJ in the following form: 

Proposition 18 Let < a < aume (where aume is defined in hS.lf^) ). then: 
(i) The semigroup T2 is uniformly exponentially contracting in Ca^ 

rxit) = sup |T2(t)M°|„ < Cexp{-K{-f,a)t)N, k>0 
uOex 
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where k{'j, a) = 7 — 



aV^ > for any ball 



X = {ue C„; 



U\a < N} 



(a) For any e > 0, the ball Ba := {m G Cq, : ImL < h e} is an absorbing set for 

bounded subsets ofCa- Here the radius a of the absorbing ball reflects the contribution of the 
free interface alone. 

Next we prove that the boundary contribution to the evolution, i.e. the operators Ti{t) 
are uniformly compact. Namely, the following proposition holds: 

Proposition 19 If a < aspace = , then for any t^ > the orbit of the ball 

^t>toTiit)X is relatively compact in Ca. 

Proof. For the version of Arzela-Ascoli theorem appropriate for it is sufficient to 
have uniform boundedness for the derivative and uniform decay of the family of functions as 
\x'\ ^00 which is faster than the decay prescribed by the weight. From the estimate (j3.18j) 
we see that the contributions from the interface decay as ex^{—{aspace + which can be 
made faster than any exp(— for a < aspace- On the other hand, the weighted estimate 
()3.25|) (cf. also Remark [TT|) demonstrate that the spatial derivative is uniformly bounded. 
Then it is easy to construct a finite e-net by choosing a finite interval beyond which the 
functions of the family are smaller than e and extending the elements of the e-net from this 
interval by zero. ■ 

The properties of the evolution operator T{t) described in the above propositions allow 
us to apply the abstract general result (see, for example, [211 Chap. 1]) that in our situation 
can be stated as follows: 

Theorem 20 The co-limit set Aa of the absorbing set Ba is a global exponential compact 
attractor for the metric space Cai Aa is the maximal attractor in Ca and it is connected. 

In order to demonstrate extra regularity of the elements of the attractor, we shall need 
a general fact concerning compact attractors (cf. [S]). 

Theorem 21 Let X be a Banach space and B G X be a ball. Let S : B ^ B be a uniformly 
continuous mapping and C G B a compact attractor for the iteration semigroup S"". Then 
on the attractor S is a mapping onto SC = C. 

The theorem holds for either continuous or discrete time. For the simplicity of presenta- 
tion we consider only the discrete case here. 

Proof. Suppose SC 7^ C, then there exists Xq G C and Xq ^ SC. In this case there exists 
a whole ball Br{xo) = {x G C : ||a; — a;o|| < t"} such that Br^Xo) H SC = 0. Indeed, since 
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the attractor is compact, its continuous image is compact and therefore the distance to xq, 
being a continuous function on a compact, attains its nonzero minimum r on SC. 

We know that S is uniformly continuous, therefore for any e there exists 6{e) so that 
ll^x — Sy\\ < e if ||x — y\\ < 6. 

Since xq € C for any e there exists n so that \\S^x — xo|| < e for some x E B, simuha- 
neously we can always select n so large that dist{S"'^^B,C) < 6. Now take y = S"'~^x and 
a G C, \\y — a\\ < 6, then \\Sa — Xo\\ < \\Sa — Sy\\ + \\Sy — Xo\\ < e + e. By selecting e < r/2 
we come to a contradiction. ■ 

As an immediate application of the above theorem we obtain the following regularity 
result: 

Theorem 22 For < a < a^j^ (where a^j^ is defined in \3.2b]) ) the semigroup on the 
attractor Aa in the space is onto; Aa consists of differentiable functions that satisfy the 
estimates ^ 

\<P\a < |0xU < Mivo, V^, a, 7) < Mivo, V^, 7) 

which yields that all Aa ^ for < a < P = a'^i^ — e. 

Proof. The theorem follows from the previous one if we note that the required uniform 
continuity follows from the well-posedness Theorem Since the mapping is onto, given 
(j) G Aa for any t there exist ip G Aa, so that (p = T{t)ip. By using estimates ()3.18p - ()3.25|) 
and taking into account exponential decay of the contribution from initial data as t ^ 00 
fnmil - (nOT|l we obtain the result. ■ 

Remark 23 Since any function in the attractor can he viewed as a result of evolution by 
the semigroup, it therefore locally satisfies the heat equation and consequently it is locally 
C°° . In addition we can show that due to the differentiability of the velocity of the interface, 
functions in the attractor are up to the interface. 

In addition to being onto, the semigroup T{t) is also one-to-one on the attractor: 

Proposition 24 The semigroup T{t) is one-to-one on the attractor Aa- 

Proof. Let T{t)ui = T(t)u2. Denote to = mi{t : (T{t)ui){x) = {T{t)u2){x)}; obviously 
to > 0. Then the difference w = T{t)ui — T{t)u2 is identically for t > to. Let ti < to 
then there exists Xo such that w{xo) 7^ 0. next we select the parabolic neighborhood (a cup) 
U{xo) = {{x,t) : ti < t < to + e, Xq — 6 — k{t — ti) < x < Xq — 5 + k(t — ti)} where e, and 6 
and k are selected in such a way that the neighborhood does not intersect the free interface. 
Since w is a solution of the heat equation in U{xo) it is real analytic in t and therefore 
should be identically in U{xo) in contradiction with the assumption w{xo) 7^ 0. ■ 

Computations for the volume evolution estimates and the Hausdorff dimension below are 
implemented in a Hilbert space. It is easy to see that the exponential decay implies the 
inclusion Ca C Ho see f|3.1U|) . Consequently, the compactness result holds for Hq as well: 
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Theorem 25 Let < a < amin- Then: 

The semigroup T2 is uniformly exponentially contracting in the H^-norm: 

^x{t) = sup ||T2(t)'U°|| < Cexp(— k(7, a)t)A^/A/a, k>0 



where ^(7, a) = 7 — — aV^ > for any ball 



X = {u e Ca] 



u\a < N}. 



For any to > the orbit of the ball Ut>taTi{t)X is relatively compact in Hq. 

The u-limit set Aa of the absorbing set Ba '■= {u E Ca '■ \u\a < \- e} is a relatively 

compact set in the HQ-metric. 

Proof. The only additional ingredient of the proof, as compared to the Ca case is 
provided by the imbedding estimate (jH.lOj) . which yields that the e-net generated for a set in 
the Ca-norm is automatically an {e / -/a)-Yiet in the i^o-norm. Since the imbedding estimate 
implies continuity of the imbedding, the above proof essentially repeats the proof of the fact 
that a continuous image of a compact set is compact. ■ 

6 Evolution of the volume elements on the attractor 

In this section we present the main result of the paper which is a proof that the Hausdorff 
dimension of the attractor is finite (for definiteness we consider ^o^i-e- « = 0). The proof 
is based on a study of evolution of the infinitesimal volume along the trajectories in the 
attractor. We demonstrate that for sufficiently large m that is defined solely by the physical 
parameters of the problem the m-dimensional volume decays exponentially. This property 
combined with the compactness suggests that the Hausdorff dimension of the attractor for 
the solutions of the free boundary problem is no larger than m. In the arguments regarding 
the Hausdorff dimension of the attractor we follow quite closely the ideas outlined in [21]. 

First we restate the problem in the coordinate frame attached to the free interface, 
X = X — s{t) as follows 



(Tildes have been omitted.) 

Let {[/(., t), \^(t)} be an orbit in the attractor. Let us consider the formal linearization 
of the problem about {U,V}: 



Ut 

g{u\x=o) 
u{x, 0) 



Uxx + v{t)ux - ju := F{u) 
v{t), [du/dx]\x=o = v{t), 
u\x). 



— 00 < X < 00, X 



(6.1) 



zt = + ZxV -^w- z{0, t)Ux/v{y{t)) ■= F'{U, V)z 



(6.2) 
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z{0,t) + u{V{t))[z,{0,t)] = 0, 



(6.3) 



z{x,0) = zo{x) (6.4) 

where i^iV) = — ((7~^)'(V^). We require I'iV) to be positive and bounded from below; I'iy) > 
uq. This condition again mimics the behavior of the Arrhenius kinetics. We have ehminated 
the velocity perturbation v{t) in the term v{t)Ux of the linearization through replacing it by 
the perturbation of the temperature z(0,t)/z/(V^(t)) that arises from the linearization of the 
kinetic boundary condition in ()6.1|1 . The linearized problem represents the first variation of 
problem (jHIH). 

It is possible to show that the linearized problem is well-posed in the following sense: 

Theorem 26 For any Zq E H there exists a unique solution z of i6.^6.4\ ) such that z G 
L2(0,T;H(t)) nC([0,T];/7) where S(t) = {/ G H\ /(O) + [/.(0)]z/(y(t)) = 0} 

Proof. This linear problem is somewhat nonstandard as it contains a nonlocal term 
(projection) z{0,t). Nonetheless it can be handled as follows. Consider first the problem 
()6.2ll6.4j) with a source, and zero initial conditions 

Wt = Wxx + WxV - 7W + J'(x, t) 

w(0,t) + [^;,(0,t)]z/(V(t)) =0, 
w(x,0) = 0, 

and let £ be its solution operator: w = £jF(x, t). Existence of unique global solutions for 
such problems is guaranteed by the general theory of linear parabolic equations. 

We regard a solution of ()6.2ll6.4j) as a superposition of an appropriate w and of W{x,t) 
which solves the homogeneous problem with the initial condition zo{x) .Then, with the non- 
local term viewed as a source, on the boundary one obtains an equation for z{0,t): 

£[-^(0, t)Uju{Vit)) + W{x, t)l=o = z{0, t) + W{0, t). (6.5) 

It is not difficult to show that the above equation is uniquely solvable as an integral equation 
with a sufficiently regular kernel. Thus, the source term is found and, consequently the 
problem ()6.2H6.4|) can be solved. ■ 

In order to estimate the evolution of the volume element we need an estimate for | |^7x| I/Zq- 
From now on, for brevity the i^o-norm will be denoted by 1 1 . 1 1 . From the imbedding estimate 
()3.1()|1 and estimate ()3.25j) we obtain the following important result: 

Lemma 27 If U E Aq, then ||f4|| < M./^/a[~ := Af, where a'^-^^ = min(^, 2^) 
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Wc arc now ready to estimate the evolution of the volume element. To this end we 
need to estimate the trace of the finite-dimensional projections of the generator of the linear 
semigroup. Let {^i, . . . , ^m} be m elements of H and let {^i, . . . , Zm} be the corresponding 
solutions of the linearized problem. Then it can be shown that the volume element spanned 
by {Ci) ■ ■ ■ ) Cm} evolves accordingly to the formula 

\zi{t) A ... A Zm{t)\ = 16 A ... A Cm)! exp / Tr [F'{U{r), ^(t)) o Qm{r)]dr, 

Jo 

where Qm{T) = Qm{T,U,V;$,i, . . . ,C,m) is the projector in H onto the space spanned by 
S(r) = {zi{t), . . . , Zm{T)} . In order to calculate the trace we need to choose a basis in S(r) 
orthogonal in the sense of H. 

Evaluation of the inner product in H gives rise to sums of integrals over the domain 
(—00, 0) U (0, 00). For brevity everywhere in the sequel we denote them by 

p pO poo 

/ f{x)dx = ( / + / )f{x)dx. (6.6) 

JR± J-00 Jo 

Let (p be an element of S(r). Consider the following inner product in H 

= -7(0,0)+ /" ci>^^(l)dx+V [ cj)^(f)dx+[(f)'{0)] [ U^(j)dx^--f+h+l2+h (6.7) 

Jr± Jr± Jr± 

We integrate Ji by parts, 

7i = 0,01° ^ + 0.01^- / cgdx 

Jr± 

It is easily seen that 

(j>^(j>dx = ]- f {(j>^ydx = 

^ Jr± 
then 

(F'0,0) = -7-[0'(O)]0(O)- / <l>ldx + [<!>' m [ UAdx (6.8) 

Jr± Jr± 

In the choice of the m-dimensional orthogonal set we shall distinguish the two possibilities: 
0(0) = (which defines an (m — 1) -dimensional subspace), and otherwise. Since the trace 
of the operator is independent of the choice of an orthonormal basis, we can choose m — 1 
basis elements satisfying the above condition. In the case 0(0) = we obtain 

(F'0,0) = -7- / 0^c^x<-7 (6.9) 
Jr± 

Note that the terms with [0'(O)] vanish since [0'(O)] = — 0(O)/z/ = in view of the boundary 
condition. 



= -[0'(O)]0(O)- / 4>ldx 

JR± 
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For the basis element with (j){0) ^ the corresponding trace component, 

0(0) 



(fijx 

R± ^ 



Ur<t>dx 



is estimated from above as follows. First we estimate the last term: 



m 



Ux(f)dx 



2a 



< TT- / Ufdx + - 



a 02(0) 



2 



(p^dx = — 
R± 2a 



U:dx + - 



a 02(0) 



2 i/2 



where a > will be chosen later on. 

Now we need the following interpolation result. By integrating from — oo to we obtain: 



^(0) 



'')xdx 



rdx < 2 



oo 



2c 



On the other hand 

J.2 



(0) 



rdx = —2 



4x<2 r(^+^)dx. 

Jo 2 2c' 



Thus 



0'(O) < 



Jr± 2 2c' Jn 



'R± 2 

In some sense, this estimate is a Sobolev trace theorem 
Finally we obtain 

-2(0) 



C02 1 

—)dx H 

R± 2 ' 2c 



< -7 + 



V 



1 



(gdx + — 



a 02(0) 



,1 a X 1 r 



l.dx + -] 
c 



2 z/2 



()2dx + — 

2a 



U^dx 



z/Q IVq I Jr± 2a J ^± 



R± 

rr9 , /I a , 1 

f/2rfx+ - + -^ - 

i/q 2i/o 2c 



If a and c are chosen so that the coefficient at the integral of 0^ is nonpositive, say if 
c = 4z/o/(2z/o + a) then 



< -7 + 



2a 



R± 



4i/n2 



Uzdx + 



for any a > 0. 

By Lemma Wj\ the norm \\Vx\\ < A/" where the bound depends only on the kinetics. 
To optimize the estimate above we choose a that gives the minimum to the expression 



)2 + 2~A/'2 considered as a function of a. This results in the estimate 



< 



fx 



. 2z/o + a 2 , 1 ^ /2z/o + 1 2 1 .r2 

-7 + mm +7r-^ 1 < ( . o ) + -A/ 



4^2 



(6.10) 
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The explicit form of /i is not important for our purposes. Thus, employing the above esti- 
mates for the trace entries ()6.9|) . ()6.1U|) we can complete the estimate for the evolution of 
the volume element: 

m 

Tr [F'(f/(r), V{r)) o Q„(r)] = 0,) < - (6.11) 

i=l 

Taking m > M = fi/j is sufficient for the trace to become negative. Note that M depends 
on z/q, 7, Vq and Vq. 

7 Differentiability of the semigroup 

To utilize the trace estimate developed in the previous section we need to demonstrate that 
the nonlinear evolution of the volume is well approximated by its linear counterpart. This 
will be ensured by the differentiability of the semigroup solving the free-interface problem 
with respect to the initial conditions, see 

For the purposes of this section we need to impose an additional condition on the kinetics 
function: we will require both g and to be twice differentiable. In applications this 
condition is definitely satisfied for all realistic kinetics. 

In Sec. 2 we cited the global existence result for the classical solutions of the free interface 
problem ()1.H) - ()1.4|1 . However, our trace estimates take place in the geometry of a Hilbert 
space. Therefore we need to introduce weak solutions by extending the existence theory to 
more general initial data that belong to a Hilbert space. The scheme of introduction of weak 
solutions is based on the following estimates: 

Proposition 28 Let U and W be two orbits (i.e., two solutions of the problem \6. 1]) ) with 
initial data Uo, Wo in the attractor: U = T{t)Uo, W = T{t)Wo. Then for any t > 

\\Uit)-Wit)\\<e'''\\Uo-Wo\\ (7.1) 

t 

I \\U{t)-Wmldr<e^' \\Uo-Wof 



where C is a uniform constant. 

Remark 29 The above proposition allows us to obtain weak solutions with initial data in the 
closure of A in the H-norm. Namely, in a standard fashion we select a Cauchy sequence of 
initial conditions in A and define the solution as the corresponding limit of smooth solutions. 
We note that it is not necessary to take initial data from A; similarly it is possible to define 
a weak solution for the initial data in the closure of a ball in Ca- 
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Proof. Let U{x,t) and W{x,t) be two solutions of the free boundary problem (in the 
frame attached to the free boundary) 

Ut = U,, + [U,{0,t)]U,-^U, g{U{0,t)) = [U,{0,t)], U{x,0) = Uo{x), 

Wt = W^, + [W,{0,t)]W^ - 7iy, g{W{0,t)) = [W,{0,t)], W{x,0) = Wo{x). 
The difference w = U — W solves the following problem 

-K(0, t)] = g{W{0, t)) - g{U{0, t)) = -{g{9)yw{0, t). 

w{x, 0) = Uo{x) - Wo{x). 

We also observe that —{giO))' < const = C while f/, W and their x-derivatives are uniformly 
bounded on the attractor. 

We multiply the equation throughout by w and integrate to obtain the following energy 
estimate for the H norm: 



ll^ll^ ~ i Wxxwdx + \Ux{Q,t)\ I Wxwdx + [wx{Q,t)\ f W^wdx — 'j f vo^dx 
= -[wx]w\q- ^\\w\\^ - \\wx\\^ + [Ux{Q,t)] / Wxwdx +[wx{S>,t)] / Wxwdx (7.2) 

JR± JR± 

We need to estimate different terms in ()7.2|) . For the first term we get on respective intervals 



K(0,t)]w(0,t)| < CM;(0,t)2 = C / {w'^)xdx<2C / \wxw\dx 



|K(0,t)]w;(0,t)| < Cw(0,t)^ j -'^^ j 1^=^^'^^ 



The sum of the above inequalities yields the estimate 

|K(0,t)]ti;(0,t)| < Cw{^,tf < C{ei \\wf + — Ww^f) (7.3) 

Next, 

\[U.{<^M I ^^^dx\ < |[?7,(0,t)]|(£2|kl|' + -|kx||') <C^i(£2|kf + -lkx||') 

JR± £2 £2 
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Also, 

K.(0,t)] [ W.wdx <C3\w{0,t)\ \\W,\\ \\w\\ 
Jr± 

< C3\w{0,t)\ Af \\w\\ < Ci{e3 \\wf + — \\w^f y/^ \\w\\ 

£3 

1 2 1 

< C'4(v^lkll + ^= W'^xW) \\w\\ = C4{^/e^\\w\\ + —— \\w^\\ \\w\\) 

< Cile^ \\w\f + — \\wr,\f) 
Collecting the estimates for different terms we get 

Ic?. .0 II ||2 II ||2 

--r\w\^ < — \\wJ\ — 7 w 

+C{ei \\wf H \\w^f) + Ci{e2 \\wf H \\w^f) + 6*5(55 \\wf H \\w^\f) 

Si £2 

/ -^11 l|2 I II ||2 

where on the last step we have chosen C/si + Cxje^ + C^/e^ < 1/2. We rewrite our last 
result as 

j^\\wf+\\wA'<C\\wf. (7.4) 

From this inequality we get first that — < C yielding, by Gronwall's inequality, 

at 

that 

< ||wo||^ exp(Ct); 
at the same time by rearranging and integrating ()7.4j) we obtain 

t t 

\wx\\^ dr < (C \\w\f ||iL'||^)(ir < llwoll^ exp(Ct). 

./ dt 



If the initial data are in fl then a similar argument yields the following estimate 
analogous to (j7.1|) : 

\\U{t)~Wmi<e^'' \\Uo-Wo\\, (7.5) 

Next we prove the differentiability in that is sufficient for the validity of the dimension 
estimate because it implies the differentiability on ^ C H^. 

Theorem 30 Let U and W be two orbits U = T{t)Uo, W = T{t)Wo, Uq, Woe H^nC^. 
Then there exists z{t) such that 



\\U{t) - W{t) - z{t)\\ < const \\Uo -WoWl 



as Wo Uo- 
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In this case the Frechet differential of T{t) at the point Uq is the mapping ^(0) = Uq — 
Wq — * z{t), where z{t) solves the linearized problem. 

Proof. The goal of the proof is to evaluate the difference between w = U — W and 
its approximation by the differential. We define z{x,t) as a solution of the free-interface 
problem linearized about the orbit U{x,t): 

zt = + [z-AO, t)]a, + [f/,(0, t)]z,, (7.6) 
z(0,t) = {g-'y{[U^{0,t)])[z^{0,t)], z{x,0) = Uo{x) - Wo{x), 

(see Theorem l26|l . For the difference y = w — z we have the following equations 

yt = Vxx + [yx{0,t)]Ux + [Ux{0,t)]y^, 
l/(0,t) = (r'y([f/.(0,t)])b.(0,t)] + ig~yie)[w,i0,t)]y2, yix,0) = 0, 

We multiply the equation throughout by y and integrate to obtain the following identity for 
the H norm: 



--II II' 
2dt 



yxxydx + [Ux{0,t)] y^ydx + [y:r{0,t)] U^ydx 
Jr± Jr± 

-[yx\y\o - hxW^ + [f/x(0,t)] / yxydx + [yx{^,t)] I U^yd 

Jr± Jr± 



(7.7) 



We need to estimate different terms in ()7.7|) 

\[yx{0,tMO,t)\<Cy\0,t) + C[w.{0,t)]'\y{0,ty < ^Cy'{0,t) + ^C[w,{0,t)]' 



Biy\0,t) + B2[w^xiO,t)Y 



< 






+^ 

2 


oo 


j yldx 


j yldx 


-oo 








2x2 



Bsisi \\wf + — \\w^\\^) 
^1 



< Bi{ei \\yf + — \\yxf) + B2.{ei \\w\\^ + — H^^;! 



2\2 



Next, 



[[/,(0,t)] / y^ydx 



< |[t/.(0, t)]| (£2 Ibf + - \\yx\?) < Cr{e2 Ibf + - hxf) 



e2 



£2 
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Also, 



[yxiO,t)] / U:,ydx 

<(C3|i/(o,t)| + i?4K(o,t)]2) / \u^y\dx<{CMM\ + B4wAo,t)]') \\y\\\\u.\\ 

< C,{e, \\yf + - \\y£Y/' \\y\\ + 55(^1 \\wf + - \\w£ + ||«;||^) \\y\\ 

< C.iV^sM + \\y.\\) \\y\\ + B,{e, \\wf + ^ \\w£) \\y\\ 

<C,{e,\\yf + ^\\y£)+B,{\M') 

where the constants C4 and include the factor \\Ux\\ ■ Collecting the estimates for different 
terms we get 

--\\yf< 

2dt - 

-\\y£ + \\yf + - + B,{e, \\wf + - \\w£)' 

El £1 

+C^ie, \\yf + - \\y£) + C,{e, \\yf + ^ \\y£) + fiedk.f) 
Now we can select ei, £2, and £3 sufficiently large so that the coefficient by is less than 



— 1/2. We collect the like terms in the above inequality to obtain 

\j^\\y\?<-\\\yA" + c\\yf + B\\w\\\ 

We rewrite our last result as 

j^\\yf + \\yA?<c\\yr + B\\w\\\ (7.8) 

from where it is clear that 

j^\\yr<C\\yf + B\\w\\\. (7.9) 
By Gronwall's inequality it yields 

t 

< B expect) J \\w\\'lexp(—Cr)dT 


<Brexp{Ct) \\wo\\t 

In the above estimate we utilized (j7.5|) . ■ 

Finally (see |21]), the estimate for the dimension of the linear volume element and differ- 
entiability of the semigroup yield the estimate for the Hausdorff dimension of the attractor: 
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Theorem 31 The Hausdorff dimension of the attractor A is no larger than 



cf. SEW. 



In conclusion it is worth mentioning that the estimate exhibits a transparent and physi- 
cally natural dependence of the dimension on the heat loss and characteristics of the kinetics 
which are the defining factors of the dynamics. We note however that numerical simulations 
|14j on ()1.H) - ()1.4|) show that the behavior without heat losses and with sufficiently low heat 
losses are qualitatively identical and exhibit the same variety of complex dynamical patterns. 



8 Acknowledgments 

The authors would like to acknowledge support in part by NSF through grants DMS-0207308 
and DMS-9704325. 



References 

[1] C.-M. Brauner, J. Hulshof and A. Lunardi, A general approach to stability in free 
boundary problems, J. Differential Equations 164 (2000), 16-48. 

[2] C.-M. Brauner and A. Lunardi, Instabilities in a two-dimensional combustion model 
with free boundary. Arch. Ration. Mech. Anal. 154 (2000), 157-182. 

[3] X. Chen and F. Reitich, Local existence and uniqueness of solutions of the Stefan 
problem with surface tension and kinetic undercooling, J. Math. Analysis Appl. 164 
(1992) 350-362. 

[4] P. Constantin, C. Foias, and R. Temam, Attractors representing turbulent flows. Mem- 
oirs of AMS, 53, no. 314 (1985) 

[5] I. D. Chueshov, Introduction to the Theory of Infinite- Dimensional Dissipative Systems, 
"Acta" Scientiflc Publishing House (2002). 

[6] S. Din, Y. Tao and H.-M. Yin, A chemical diffusion process with reaction taking place 
at free boundary, Canad. Appl. Math. Quarterly 5 (1997) 49-74. 

[7] M. Frankel, M. Qu and V. Roytburd, On a free interface problem modeling solid com- 
bustion and rapid solidiflcation in inflnite medium. World Scientific Series in Applicable 
Analysis, Vol. 4- Dynamical Systems and Applications, R. P. Agarwal ed.. World Sci- 
entiflc (1995), 263-278. 



31 



[8] M. Frankel, V. Roytburd, and G. Sivashinsky, A sequence of period doubling and chaotic 
pulsations in a free boundary problem modeling thermal instabilities, SIAM J. Appl. 
Math. 54 (1994), 1101-1112. 

[9] M. Frankel V. Roytburd G. Sivashinsky, Complex dynamics generated by a sharp inter- 
face model of self-propagating high-temperature synthesis. Combust. Theory Modelling 
, 2 (1998), 479-496. 

[10] M. Frankel, G. Kovacic, V. Roytburd, and I. Timofeyev, Finite-dimensional dynamical 
system modeling thermal instabilities, Physica D 137 (2000), 295-315. 

[11] M. Frankel and V. Roytburd, Finite-dimensional attractors for a free-boundary problem 
with a kinetic condition, Appl. Math. Lett. 15 (2002), 83-87. 

[12] M. Frankel and V. Roytburd, Compact attractors for a Stefan problem with kinetics, 
EJDE (2002), 1-27. 

[13] M. Frankel and V. Roytburd, Finite-dimensional attractor for a one-phase Stefan prob- 
lem with kinetics, J. Dynam. Differential Equations 15 (2003), 87-106. 

[14] M. Frankel and V. Roytburd, On attractors for a sharp interface model of exothermic 
phase transitions. Advances Math Sciences Applications 14 (2004), 25-40. 

[15] A. Friedman, Partial Differential Equations of Parabolic Type., Prentice-Hall, Englewood 
Cliffs, N.J. (1964). 

[16] S. M. Gol'berg and M. I. Tribelskii, On laser induced evaporation of nonlinear absorbing 
media, Zh. Tekh. Fiz. (Sov. Phys.-J. Tech. Phys.) 55, 848-857 (1985). 

[17] J. S. Langer, Lectures in the theory of pattern formation, in: Chance and Matter, J. 
Souletie, J. Vannimenus and R. Stora, eds., Elsevier Science Publishers (1987). 

[18] L. Lorenzi, A free boundary problem stemmed from combustion theory. Part I: Exis- 
tence, uniqueness and regularity results, J. Math. Anal. Appl. 274 (2002) 505-535. 

[19] S. Luckhaus, Solutions for the two-phase Stefan problem with the Gibbs-Thomson law 
for the melting temperature. Euro. J. of Appl. Math. 1 (1990) 101-111. 

[20] B.J. Matkowsky and G. I. Sivashinsky, Propagation of a pulsating reaction front in solid 
fuel combustion, SIAM J. Appl. Math. 35 (1978), 230-255. 

[21] Z. A. Munir and U. Anselmi-Tamburini, Self-propagating exothermic reactions: the 
synthesis of high -temperature materials by combustion. Mater. Sci. Rep. 3 (1989), 
277-365. 

[22] E. V. Radkevich, The Gibbs-Thomson correction and conditions for the existence of a 
classical solution of the modified Stefan problem, Soviet Math. Dokl. 41 (1991) 274-278. 



32 



[23] K. G. Shkadinsky, B. I. Khaikin and A. G. Merzhanov, Propagation of a Pulsating 
Exothermic Reaction Front in the Condensed Phase, Combust. Expl. Shock Waves 7 
(1971), 15-22. 

[24] R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physics, 
Springer Verlag: New York (1988). 

[25] Van Saarloos, W. and Weeks, J. , Surface undulations in explosive crystallization: a 
nonlinear analysis of a thermal instability, Physica D 12 (1984), 279-294. 

[26] A. Varma, A. S. Rogachev, A. S. Mukasyan and S. Hwang, Combustion synthesis of 
advanced materials. Adv. Chem. Eng. 24 (1998), 79-226. 

[27] A. Varma, Form from fire, Sci. American Aug. (2000), 58-61. 

[28] la. B. Zeldovich, G. I. Barenblatt, V. B. Librovich and G. M. Makhviladze, The Mathe- 
matical Theory of Combustion and Explosions, Consultants Bureau, Washington, D.C., 
1985. 



33 



